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TECHNICAL NOTE 2853 


A STUDY OF THE APPLICATION OF POWER-SPECTRAL METHODS: 
OF GENERALIZED HARMONIC ANALYSIS ТО. 
GUST LOADS ON AIRPLANES 


By Harry Press and Bernard Mazelsky 
SUMMARY 


The applicability of some results from the theory of generalized 
harmonic analysis (or power-spectral analysis) to the analysis of gust 
loads on airplanes in continuous rough air is examined. The general 
relations for linear systems between power spectrums of а random input 
disturbance and an output response are used to relate the spectrum of 
airplane load in rough air to the spectrum of atmospheric gust velocity. 
The power spectrum of loads is shown to provide a measure of the load 
intensity in terms of the standard deviation (root mean square) of the 
loed distribution for an airplane in flight through continuous rough air. 
For the case of а load output having a normal distribution, which appears 
from experimental evidence to apply to homogeneous rough air, the standard 
deviation is shown to describe the probability distribution of loads or 
the proportion of total time that the load has given values. Thus, for 
an airplane in fligbt through homogeneous rough &ir, the probability 
distribution of loads may be determined from а power-spectral analysis. 


In order to illustrate the application of power-spectral analysis 
to gust-load analysis and to obtain an insight into the relations between 
loads and airplane gust-response characteristics, two selected series of 
calculations are presented. In the first series, the standard deviations 
of lcads in continuous rough air described by an assumed power spectrum 
are calculated for systematic variations in the frequency ani damping 
characteristics of the airplane response to a step-gust input. The 
results obtained indicate that the loads in rough air are particularly 
sensitive to variations in the damping characteristics of the oscillatory 
response to а step gust and largely independent of variations in the 
frequency. In the second application, the standard deviation of loads 
is calculated for selected variations of each of several airplane geo- 
metric and aerodynamic parameters of an idealized and stable transport- 
type airplane. The standard deviations obtained are compared with results 
derived by conventional techniques of using the calculated peak response 
to an idealized and representative discrete gust. The results indicate 
that for stable configurations both methods of analysis yield results 
that are consistent to a first approximation. 
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INTRODUCTION 


The study of gust loads on airplanes 3в a twofold problem requiring 
the adequate representation of the characteristics of atmospheric turbu- 
lence and the determination of the airplane response (loads or motions) 
in rough air. These problems have teen recognized since the inception 
of gust-load research but because of the difficulties involved only 
limited approaches to the problems appeared practical. The methods that 
have been used are described and discussed in detail in reference 1. In 
general, the approach has been to use simplified airplane-response theory 
for the determinaticn cf the characteristics of discrete gusts from air- 
plane measurements of load. The "guste" derived on this basis are then 
used to calculate loads on other airplanes. Although these procédures 
appear reasonable for transferring loads to similar airplanes, ав indi- 
cated in the reference, they are of questionable value for airplanes of 
widely different characteristics (such as, configurations and stability 
characteristics). These limitations have, however, not proved serious 
in the past since the transport airplanes which were primarily affected 
by gust standards appeared, in general, to follow conventional design. 
Available data indicate that new transport airplanes experienced gust 


loads which were, in general, compatible with those predicted from past 
work. 


Trends in aeronautics. toward higher speed and the deveiopment of 
missiles have served to introduce a widening range of unusual configura- 
tions and aircraft stability characteristics. Furthermore, the gust- 
load design requirements, which formerly were of concern for transport 
and bomber airplanes only, appear to have become important for other 
aircraft as well. As a consequence of these developments, the need for 
more generally applicable techniques, both for the measurements of the 
characteristics of atmospheric gusts and for the calculation of the gust 
loads on new airplanes, has become more urgent. 


Developments in the theory of generalized harmonic analysis 
(ref. 2) appear adaptable for extending the analysis of gust loads beyond 
the discrete-gust case to the case of continuous turbulence. Techniques 
from generalized harmonic analysis involving the concept of power spec- 
tral density have been used for many years in diverse fields, such as 
іп the study of random-noise problems in communications and in the study 
of small-scale turbuience of vind tunnels. “Тһе concept of the power of 
a random disturbance, which is fundamental to the present study, is 
defined by analogy to electrical power to be the time average of the 
square of the disturbance. The portion of the power arising from com- 
ponents having harmonic frequencies between w and 0 + do is denoted 
as the power spectral density. For linear systems, the power-spectral- 
density functions of & random input disturbance and an output response 
are related through the frequency-response characteristics of the system. 
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Attractive features of spectral analysis for the study of gust 
loads are the possibilities that: 


(1) Continuous turbulence can be described in analytic form by 
а power spectrum rather than by discrete gusts. 


(2) The load response of airplanes to continuous rough air can be 
evaluated. 


(3) The desirable response characteristics of an airplane for 
minimizing gust effects in continuous rough air will become amenable to 
analysis. 

In view of the attractive features of power-spect:al-density 

- methods of analysis, an investigation of the applicability and implica- 

tions of these techniques to gust-load analysis was undertaken and the. 
results obtained are reported herein. In this paper, the concepts and. 
relations of generalized harmonic analysis are defined and their method” 
of application to the gust-load problem is indicated. The applicability 
of the normal probability distribution for the representation of the и 
probability distribution of loads in continuous rougk air is considered. 
For the case of a normal distribution of loads, the standard deviation 
(root mean square) of the load history defines the probability distribu- 
tion of loads. It is indicated that the standard deviation may be deter- 
mined from load power spectrum; thus the resuits of a power-spectral 

T analysis permit the determination of the probability distribution of 
loads. Finally, in order to illustrate the application of power-spectral 
analysis to gust-load calculations and to obtain an insight into the 
relation between loads in continuous rough air and airplane character- 
istics, two applications are presented. Іп both applications, the power 
spectrum of atmospheric turbulence obtained from flight measurements 
(ref. 2) is used to represent the turbulence input. The first applica- 
tion is intended to represent the effects on gust loads of variations in 
airplane dynamic longitudinal stability. The second application is 
intended to illustrate the effects on gust loads of variations of some 
geometric and aerodynamic parameters of ап idealized transport airplane. 
Ihe indicated variations in load intensity are compared with those derived 

= by conventional techniques of using the peak-load response to an idealized 
discrete gust. : 
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SYMBOLS 


n А( ) - response to unit step disturbance 
а+- а slope of tail lift curve per radian 
\ i : 
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ct 


T( ) 


slope of wing lift curve per radian 
parameters of equation (31) 
airplane wing chord, ft 

airplane tail chord, ft 
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pitching moment of inertia about center of gravity, slug-ft= 


reduced frequency, ос/У, radians/chord 


- 


horizontal distance from center of gravity of airplane to 


wing aerodynamic center, ft 


horizontal distance from center of gravity of airplane to 


tail aerodynamic center, ft 
number, of observations 
acceleration increment, g 


probability 


- 


cumulative probability distribution of z, equation (28) 


distance, chords 

wing area, sq ft | 
tail area, ва ft 

time, sec 


arbitrary value of t, вес 


frequency-response function with respect to argument, 


о, 0, or У 
vertical gust velocity, ft/sec | 
airplane true airspeed, ft/sec | 
airplane weight, 10 | 
response to unit impulse disturbance 
distance, ft | 


arbitrary value of x, ft 
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de/da 


Subscripts: 


max 


basic 


arbitrary random variable 

probability density function of variable y, equation (23) 
arbitrary function of arguments 5, x, and s 

average power of y(t), equation (4) 

standardized variable, An/o 

coefficient of skewness, equation (30b) 

coefficient of kurtosis, equation (30c) 

standard deviation of specified random variable, equation (2b) 
mass air density, sluge/cu ft 

time displacements, sec i 

autocorrelation function, equation (9) 


power-spectral-density function of an arbitrary disturbance 
with respect to о or Я i 


power-spectral-density function of a designated input 
power-spectral-density function of a designated output 


normal distribution with mean of О and ptandard deviation 
equal to 1, equation (27) 


circular frequency, radians/sec . 


reduced frequency, o/V, radians/ft 


downwash factor 


maximum ова response to а discrete gust 


. basic airplane configuration 
A bar over a symbol designates the average value of the quantity. 


In this analysis, the use of several independent variables +t, x, 


and в for an arbitrary disturbance у( ) and their associated fre- 
quency arguments о, 0, and V has been found necessary. In order to 
designate that the several functions y( ), Ф( ), апа T() depend 
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upon their arguments а circumflex ~ and а tilde ~ have been used . 
over the appropriate sets of functions in accordance with the scheme 


shown in the following table: 
Variable | 
Function 


Frequency argument 


Disturbance 


Power-spectral-density function 


Frequency-response function 


Impulse response 


Step response 


ANALYSIS 


Basic Concepts and Relations of Power-Spectral Analysis 


The theory of generalized harmonic analysis is an outgrowth and 
generalization of harmonic analysis and is largely the work of Norbert 
Wiener (ref. 3). Accounts of the theory snd method of application are 

. also given in references 2, 4, and 5. Іп order to orient the reader, a 
brief account of the background, the basic concepts, and the relations 
4 is presented. 


Тһе theory of harmonic analysis indicates that an arbitrary periodic 
function can be.represented by a Fourier series in the following manner: 
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where 


ta =È | PO) cos nat at 
2 pt | 
Ba =F | F(t) sin not dt 


and Т is equal to 2л/® and is the period of the function. In order 
to apply this technique to nonperiodic phenomena, the limit of equa- 
tion (1) as T goes to infinity must be considered. For this case of 
a nonperiodic function, the Fourier series takes the form of the Fourier 
integral and is given by 


F(t) = 2. a еї Г ^ F()e ат (2) 


If the second integral on the right is denoted as G(a) then equa- 
tion (2) can be written in the reciprocal form 


F(t) = EI G (o) e “д © (3) 
where 


С (a) -Г F(1)e Уа 


со 


The Fourier transform pair thus provides reciprocal relations between 
the time function F(t) and its frequency representation G(m). The 
quantity G(w) dm gives the contribution of those harmonic components 
of F(t) whose frequencies lie between œ and о + до, | 


A necessary condition for the application of equation (3) is that 
the integrals involved be convergent. This condition acts as a severe 
limitation on the applicability of the Fourier integral relations. Іп 
many problems, such as neise in an electric circuit or turbulence encoun- 
tered by an airplane in flight, the disturbance is nonperiodic, persists 
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for a relatively long period of time, and shows no tendency of dying 

out. Іп these cases the required integrals do not converge and, as ~ 
а consequence, the frequency representation in terms о? equations (2) | 

and (3) is not possible directly. 


In order to develop a frequency representation which would be 
. applicable to continuing disturbances, the theory of -random processes 

makes use of the concept of a stationary random process. The character- 
istics of a stationary random process are described in detail in refer- 
ence №. Essentially, the assumption is that the underlying mechanism 
which gives rise to the disturbance does not change in time and that a 
statistical equilibrium exists. Thus, the statistical characteristics 
of the distribution are invariant with time and statistical prediction 
becomes possible. For the case se of а stationary random function of 


time y(t), the mean square у2(%) is defined by 


Й | | | 
Е _ -y2(t) = Qum di pas | (5) | 


— 


Ы ын а ЧИ Ри па ТУ 
* 


- 


The mean square will usually exist and represent a measure of disturb- 


ance intensity. Since у2(%) ів а quadratic function of y(t), it has | 
been termed the "average power" of y(t) in analogy to electrical power 

which is proportional to the square of the current. The function y(t) 

is considered to be composed of an infinite number of sinusoidal сош- 

ponents with circular frequencies ш, between О and о, The portion 


of у2(+) arising from components having frequencies between © and 

о + do is denoted herein by 6(m) дю. Тһе function Ф(о) has beencalled 
the power-spectral-density function in the literature. From this defini- 
tion, (œ) has the property that | 
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| j ; - 
| The power-spectral-density function of a галаст variable y(t) is 
generally defined in the following manner (see for EM ref. 4) 
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where the notation | | indicates the modulus of the. сопрјех quantity. 
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A ЧОН 


If Ф(о) is defined in this manner, it has the property that its 

integral over the limits О to » 18 equal to the power т. That 
this expression for (w) is consistent with the preceding discussion | 
is seen to be plausible from the fact that, for each frequency, Ф(о) | 
is proportional to the square of the amplitude of the component of y(t) | 
at that frequency. Thus, Ф(о) dm is a measure of the contribution of 

that frequency to y(t). It should be mentioned that the definition | 
equation (б) differs in minor. detail from that used in reference 2 but | 
agrees substantially with those used in references 4 and 5. . 


А significant and useful relation for linear systems exists between . i 
the power-spectral-density function %;(m) of a random input disturb- - 
ance and the power-spectral-density function ® (œ) о? an авеосіафей 1 
output through the system frequency-response function. T(io). The жамы 
system frequency-response function (or admittance) Т(1о) is defined } 


such that Т(їш@)е1®© ‘ів the system response to the sinusoidal input eiat, | | 
In these terms, the relation between the power-spectral-density functions 
18 given by ; 


dla) = dla) [TC] 0 7. 


For a given linear system, the function T(im) may ђе conveniently. 
obtained from < „nur the unit impulse reaponse ог the unit step response, 
respectively, by means of the following relations: 


. T(1«) -f W(t)e^ 19 tag А | 


(8) - 
T(im) = to [alte at | | | 


where W(t) is the response to а unit impulse and A(t) is the response 
to a unit step. i - | 


Equation (6) may be used to evaluate the power-spectral-density 
- function from observed data. However, in practice, the power-spectral- 
density function may be determined more conveniently апа less tediously 


| 
| 


10 | NACA TN 2853 


~ 


by using a related function, the autocorrelation function R(t), 
defined by А: 


R(t) = lim = " y(t) y(t +7) at (9) 
Т —> со 0 


The autocorrelation function has the symmetrical property R(t) = К(-т). 
and is reciprocally related to the power-spectral-density function by 
the Fourier cosine transformation in the following manner: 


В(т) d lw) сов ‘ar do 


(10) 


со 
Е òla) = a R(t) cos wr at 
0 : 


Reference №, for example, shows that this definition of the ромег- 
spectral-density function is consistent with the preceding definition 
of equation (6). | 


Some Forms of Power-Spectral Relations for Gust Applications 


In this section, the method of application to gust loads on air- 
planes of the aforementioned concepts and relations is considered. Some 


available information is presented on the frequency-response function  T(io) 


for the gust-load condition and on the power-spectral-density functions 


for atmospheric turbulence. Finally, for convenience in -gust-load appli- 


cations of power-spectral methods, the functions and relations are pre- 
sented in terms of distances rather than time. Use is made of the dis- 
tances x in feet and s in chords and their associated frequency 
argurents Я and v. These changes of variables are subsequently shown 
to be particularly appropriate in the gust case because in these terms 
the power spectrums are independent of airplane forward speed. i 


If atmospheric turbulence can ђе considered a stationary random 
process, then the basic requirement for the application of the foregoing 
power-spectral concepts and relations is satisfied. Evidence that this 
assumption is plausible under some conditions exists and is discussed 
subsequently. On the basis of this assumption, the turbulent-vertical- 
velocity cistribution along а line in space сап at an instant in time 
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be considered to represent а stationary random function of space cone 
sisting of an infinite number of harmonics of various frequencies or 
wave lengths. For this condition of spacial wave lengths, a natural 
unit for the associated frequencies would appear to be radians per unit 
distance or radians per foot. The description of the power-spectral- | 
density function of atmospheric turbulence must thus be given basically 
in terms of units such as radians per foot. However, for an airplane in 
flight through rough air, consideration of the airplene as penetrating 
the gusts and experiencing the associated loads in terms of time is fre- 
quently convenient. Thus, the power-spectral-density functions of gust 
velocity and loads experienced by the airplane may be considered expressi- 
ble in terms of the frequency argument œ in radians per second used in 
the preceding section. It is, therefore, permissible to express the 
relation between the power spectrums of gust velocity and loads or normal 
acceleration in terms of equation (7) as 


dla) = $; (в) |T]? 


where 0 (œ) is the power-spectral-density function of airplane normal 
acceleration, Фі(0) is the power-spectral-denaity function of gust 


velocity experienced by the airplane, and T(io) is the airplane normal- . 


acceleration response function for a sinusoidal gust velocity input. 


The present study is primarily concerned with vertical gust velocity 
inputs and acceleration increment An (load factor) outputs. Тһе func- 


tion Фі(о) will have the dimensions of (ft/sec) 2/radian/secs |T( 10)| 2 
will be given in (g/tt/sec)?. Consequently, o(a) will have the 


dimensions of g°/radian/sec. In this form, the power-spectral-density 
functions, heving the dimensions of power/radien/sec, refer to a par- 
ticular airspeed. Before considering the representation of these func- 
tions іп a form independent of airplane speed, some гешэгкв on the 
determination of the frequency-response function for the gust-load con- 
dition and on available information concerning the power-spectral-density 
function for atmospheric turbulence appear appropriate. 


Frequency-response function.- Experimentsel methods for determining 
the frequency-response function of an airplane for а gust velocity dis- 
turbance have unfortunately not yet been developed. The frequency- 
response function for airplanes for a gust disturbance can, however, be 
estimated by theoretical methods by solution of the airplane equations 
of motion for a gust disturbance. Methods for the determination of the 
frequency-response function for a linear system are described in chapter 2 
of reference 4, for example. The calculation of the frequency-response 
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function may sometimes be more conveniently performed by first deter- 
mining the airplane response to a step gust by methods such as those 
described in reference 6. The response of а step gust may then be 
used in equation (8) in order to determine the frequency-response 
function for a continuous sinusoidal gust input, 


3 Power spectrum of atmospheric turbulence.- The power-spectral- 
density functions of atmospheric gust velocity have been studied from 
airplane flight measurements (ref. 2). Measurements of pitching velocity 
of a B-25 airplane in flight through rough air were used to determine 

the power-spectral-density function of atmospheric turbulence for four 
weather conditions. Clementson determined the reduced autocorrelation 
function R(T)/R(O) of pitching velocity and obtained the output power- 
spectral-~density function by taking its Fourier transform. Тһе appro- 
priate frequency-response function as determined from simplified theoreti- 
cal calculations was then used in equation (7) to obtain the input spec- 
trum, The normalization of the autocorrelation function R(t) by dividing 
through by R(O) is arbitrary and was made on the basis of yielding a 
pitching-velocity power of 1 (Gegree/second)? for tbe airplane used in 
that investigation. The normalized power-spectral-density functions of 
atmospheric turbulence derived in this manner did not vary appreciably 
between weather conditions; although, as might be expected, the value 

of R(O) for the pitching velocity output did vary appreciably, and 
reflected variations in the average power of turbulence with weather 
condition. Оп the basis of the results obtained in reference 2, the con- 
clusion was reached that "atmospheric turbulence is a stationary random 
process that can be statistically described by a single reduced power- 
spectral-density curve." Although additional tests under a wide variety 
of atmospheric conditions and the use of other airplanes are needed to 
verify this conclusion, the spectrum obtained appears to be representa- 
tive at least of the conditions covered by the tests. .The results 


.obtained in reference 2 for the power-spectral-density function of 


atmospheric turbulence thus provide a turbulence input and are used in 
subsequent applications in the present study. 


The average reduced power-spectral-density function obtained in 
reference 2 for an airplane true flight speed of 300 feet ‘рег second is 
shown in figure 1. The results shown are corrected for two errors made 
in reference 2 and pointed out in reference 7. In addition, the results 
of reference 2 have been divided by 2x т order to conform to the 
definition of the power-spectral-density function used herein, equa- 


tion (6). The power-spectral-density function shown in figure 1 has 
dimensions of (ft/sec) /radians/sec and applies to an airplane airspeed 
of 300 feet per second. 


Changes іп the frequency argument.- As pointed out previously, Ф(а) 
has the dimension of a power/radian/sec and thus depends upon the airplane 
forward speed. Іп order to express the power-spectral-density functions 
in terms independent of airplane flight speed, the change of variables 
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= Q | 
(11) i | 
x = Vt J 
is introduced. The variable 0 is a reduced frequency in radians per 
foot. The variable x ів the airplane flight distance in feet, In - 
terms of these variables, the average power of a disturbance f(x) is 
given by = B › 


a р қ * ` 2 = 
92(х) = lim ы ¥(x) | ах (12) 
у - 0 Е ) , - E 


(The use of the circumflex is explained in the section "Symbols.") 


If the power-spectral-density function Ф(0) 1s given by | 


а x И > = 
Е | $(8) = lim x Ј, уже axl - (13) 


then 6(2) gives the average power of the disturbance arising from com- 
ponents having frequencies between 2 and 0 + 40. ‘Inasmuch ав 


-— = 


| ва) КС и E ; (15) E 
о о _ i 


= 


it follows from equation (11) that 
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о(а) = Уо) = volva) _ 4 _ (15) 
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; As pointed out in reference 7, the power-spectral-density function 
| Е for atmospheric turbulence representing $,(9) ів incorrectly given in 


: reference 2. The power-spectral-density function $4 (0) which corresponds 
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to the function shown in figure 1, was obtaincd by equation (15) and is 
shown in figure 2. The resulting power spectrum of atmospheric turbu- 
lence is independent of airplane speed and in a limited sense may repre- 
sent a universal reduced power-spectral-density function of atmospheric 
turbulence., | 


The relation between the frequency-response functions for the 


-sinusoidal input eiat апа е10Х {в now considered. Since for a 


given airspeed the gusts represented by elOt апа еї9х are, from 
equations (11), the same gusts, the load responses when expressed in 
terms of time t and distance x will only involve a change of scale. 


If Т(і0)е10Х ің the airplane load response to the gust e4X, then 
the frequency-response functions are.simply related in the following 
manners 


(19) = T(io) = T) "^ (16) 


The frequency-response function T( 10) may be obtained from relations 
similar to those of equation (8) and given by 


$(10) = L W(x) e7 12% gx 


i (17) 


^ mx 10 
T(in) = ia || А(хје “ах 
0 


where W(x) and A(x) are the response functions of the unit impulse 
апа unit step disturbances, respectively, expressed as a function of 
flight distance. These response functions are related to W(t) апа 
A(t), respectively, by the relations 

W(x) = 5 w(t) 

A(x) = A(t) 


In terms of the variable 0, the input-output relation can, from 
equations (15) and (16), be expressed by the relation 


Фо(2) = 35 (2) |т( 10)/ | (18) 


b 

р 

+” — -A ттан нан саша ны нене MALA карарын пола 
LI 


7^9 k "oai o 


т <<“ -w -ь..---- 


A ~ тыб а — „= л ee —. 
[4 


7. -чь У и 
Г] 


NACA TN 2853 | | 15 


"The average power for & given input and output in terms of the fre- 
quency N may from equations (1h) апа (18) be obtained ty the relations 


Pea) = Јр 8.48) an = "або an a» 


Another set of variables is useful in the applications of power- 
spectral analysis to gust loads. This set involved the reduced fre- 
quency У in radians per chord and the nondimensional distance в іп 
chords. Because of the need for consideration of unsteady aerodynamics = 
in gust-loads analysis, calculated responses to step gust velocity inputs 
may frequently be given in terms of the nondimensional distance в and 
the frequency-response function in terms of the reduced frequency v. 
These variables V апа в are related to о and 0 and to t^ end x 
in the following manner; 


‘ 
1 праны ee ^de. Rte qn di n 


(20) 
v = 5 = са 


If Т(іу)е!У8 ів the acceleration response to the gust е1У8 then 


Т(іу) may also be obtained from the impulse or step response by the 
relations 


Tv) = | " #(в)е-178ав 
(21) 


(1) = iv Г А(з)е-1УЗ ав 
0 


where Я(5) апа A(s) аге the responses to the impulse and step gust 
inputs in terms of the variable в. These response functions are related 
to W(t) and A(t), respectively, by the relations 


W(s) = 7 w(t) 


A(s) = A(t) 
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The frequency-response function (37) is related to the func- 
tions Т(іш) and T(iN) in the following manner: 


Т(10) = #(1с0) = т(ауп) (22) 


and involves only a change of scale of the frequency axis. These геја- 
tions are used in the subsequent applications. 


- 


Relation of Power-Spectral-Density Function to Applied Gust Loads 


In order to apply the foregoing methods of analysis to the study of 
gust loads, the power spectral density of loads must be related to the 
intensities of the actual loads. The power-spectral~density function, 
which provides a measure of the average power arising from components 
at various frequencies, does not directly reflect the load intensities 
since the actual load at a given time represents the combined output at 
the various frequencies. Thus, it is desirable to relate the power 
spectrum of loads to specific quantities of concern in load studies, 
such аз the proportion of total time at a given load intensity (proba- 
bility distribution of load intensity), the number and intensity of peak 
loads, and other such particular quantities that may be of interest for 
structural analysis. Іп the present report, only the relation of the 
power-spectral-density function to the probability distribution of load 


is considered and more specifically the significance of the normal dis- 
tribution for loads is considered. 


Probability distribution of output.- When а linear system 18 
exposed to an input varying in a random manner with time, the probability : 
distribution of the system output y сап frequently be represented by a 
normal probability density distribution defined by the relation 


(23) 


where У and с are the mean and standard deviation and are defined by 


T 
Y= lim i] yat (оца) - 
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= T 4 1/2 
с = | lim iJ. (y - y) at (24b) 
T —» го 0 


Investigations of communication problems associated with noise, which 

has many obvious similarities to turbulence, have shown that normal dis- 
tributions are frequently encountered. Rice in reference 8 has, for 
example, shown that for a linear system the shot effect in a vacuum 

tube gives rise to a noise current which has a normal distribution of cur- 
rent intensity. Investigation of fluid turbulence frequently yields 
normal distributions of velocity fluctuations. These results appear to 

be explained by the central limit theorem of probability (ref. 9) which 
states that, under general conditions, the distribution of the sum of a 
large number of random variables tends toward a normal distribution. 


The theoretical derivation of the probability distribution of loads 
for an airplane in flight through rough air is mathematically difficult 
_ апа appears to involve assumptions regarding the nature of turbulence 
that are, for the present, questionable. However, if the load time history 
in continuous rough air is considered to be a stationary random function, 
the time history сап be considered to be made up of the sum of a large ` 
number of harmonic components with random phases, ІТ the associated 
power-spectral-density function is relatively uniform (having no sharp 
peaks), then the load at any given time is the sum of a large number of 
random variables of roughly the same order of magnitude. These condi- 
tions meet the principal requirements for applicability of the central 
limit theorem and, consequently, it follows from tbis theorem that the 


probability distribution of load intensity may tend toward a normal 
ВНЕ ОВ 


Significance of the normality of the distribution of loads.- If the 
distribution of airplane acceleration increment in continuous rough air 


is normal with а zero mean value, ihe probability density distribution 
of the acceleration increment ^n is completely described by Mic standard 
deviation and from equation (23) is given by 


flan) = —-e | | (25) 
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uation (25) is termed the probability density function of An, and 
| Tian) а(Ап) can be considered to represent the proportion of total 
E tius (or total flight distance) that An has a value between Ап 
апа Ап + d(An). 
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« 
. 


. In order to examine the properties of equation (25), it is con- | 
J venient to corsider the distribution of the тате | 


z=% В (26) _ 


Where z is the so-called standardized variable and has the probability 
Ж density distribution Е 


| -22/2 и _ 
tz) == e"s (D - 
| | Я | х _ | | 


- , oe ' 1 


The function У(2) is the normal distribution with a mean of zero and а | 
standard deviation equal to 1 and is the error function commonly tabu- С | 
Ж lated. The probability that a random value of z will exceed а given | 

value is given by the integral of equation (27) as follows: © 


-- 


ға) = f " v(z) аз | — (eB) | 
2 сы 


~ 


Equation (28) defines the cumulative probability distribution of the 
reduced variable 2 and in the case of time-history data may be con- 
sidered to represent the proportion of total time that the value of 2 
с exceeds a given value. For a fixed value of probability Р, 2 

is fixed and can be obtained from tabulations of the integral of the 
error function. Thus, for example, for Р = 0.00275, 2 = 2. 


Equetion (28) indicates that the probability of exceeding a-given 
value of 2 is a function of only the given value. Conversely for a 
given value of probability, the largest value of 2 exceeded is also Е 
fixed and depends only upon the probability. Making the substitution 


2 = ба. into equation (28) yields the result that the probability of - 
exceeding а given value of An/o is likewise a function only of the 
value of An/o and is given by 2035 The. converse also applies that, Н 
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for a given value of probability, the largest value of An/o exceeded 
is fixed and depends upon the probability. Consequently, the largest 
value of An exceeded with a given probability depends upon the value 
of с and is given by oz. The largest value of An .exceeded with a | 
given probability is thus seen to be directly proportional to the 
standard deviation. For example, when the standard deviation is 
doubled, the largest value of An exceeded with a given probability 
is also doubled. This direct relationship between the value of the 
standard deviation and the load probability distribution makes the 
standard deviation an important and significant measure of the load 
experience for the case of gust loads having a normal distribution. 
Because of this direct relation between the standard deviation and the 
load probability distribution, the standard deviation will be useful 
as a measure of load intensity in the present study. 


The foregoing discussion has served to establish the plausibility of р 
gust loads having a токтај distribution under some conditions and the sig- 
nificance of the standard deviation of loads for the case of a normal dis- 
tribution. The standard deviation of loads can be derived from a power- | | 
spectral analysis in the following manner: From its definition (eq. (24b)), 
the standard deviation of load increment output having a zero mean value 7 
is the square root of the average power. Thus, the росчег-вресіга1-2епві%у 4 
function of loads апа the standard deviation of the probability distri- 
bution of loads are related. The standard deviation o may thus be 
obtained directly from the power-spectral-density functions by the 


relation у р 


M 5 6 (m) am = [* 209) an |. (eg) | 


This relation between the probability distribution and power-spectral- 
density function for the case of a normally distributed output ties the 
power spectrum to a basic characteristic of the load history and is 
thus of importance for applications to gust-load analysis. 


APPLICATION TO GUST-LOAD PROBLEMS 


In view of the simple relations between the loads and the power- 
Spectral-density functions for the case of normal distributions of 
loads, the determination of the normality of load distributions appears 
to be an important problem in the application of power-spectral analysis 
to gust-load problems; therefore, some experimental gust-load distribu- 
tions are examined for normality. Two applications of power-spectral 


С 5у 5 А | E 
Cun». = - Iq. 


methods to the calculation of gust loads are also given in this section 
in order to obtain some insight into the relations between airplane 
characteristics and loads in continuous rough air. . In the first series 
of calculations, the input-output relations are used to calculate the 
standard deviation of loads for a selected series of idealized responses 
‘to a step gust. The second series of calculations is made to determine 
the variations in the standard deviation of load for individual varia- 

° tions of certain geometric and aerodynamic parameters of an idealized 
transport airplane. The indicated variations in load intensity in this 
application are compared with those derived by conventional techniques 
of using the peak load response to an idealized discrete gust. 
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и Observed Distributions of Loads | : 


In order to determine whether frequency distributions of load are 
actually normal distributions, recourse мав made to some available 
experimental gust-load time-history data. Time-history records of the . 
normal acceleration for two airplanes of thé same type (differing only ` | 
slightly in center-of-gravity position) in side-by-side flight through ; 

-continuous rough air at low altitudes above generally flat terrain were | 
available from a recent investigation. The test conditions and summaries ' 

of the data are presented in reference 10. The normal-acceleration 
time-history records for about a 2-minute section of one run at 450 miles 
per hour were evaluated in detail by taking readings at intervals of 
1/20 second, roughly one reading for each 4 chords of airplane travel. 
The frequency distributions for the two airplanes are summarized in 
table I. (Additional data of the same type were also examined but do 
net appear to warrent reporting in detail at this time.) In determining 
these distributions, the 1 g level was assumed to be at the mean value 

of the disturbance. (This assumption is frequently used in gust evalua- 
tions in view of the difficulty in exactly determining the 1 g line for 
the flight condition.) The other primary characteristics of the observed 
frequency distribution, standard deviation с, coefficient of skew- 

ness аҙ, and coefficient of kurtosis œ, are also given in table І 


$ 
| 
where the statistical characteristics of the distribution were determined | 
by the following relations: 
1 


N 


5 1/2 | 
o pe = 592 | (30a) 
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where N is the number of observations. For a normal distribu- 
tion аз = О and а, = 3. The distributions of table I are shown ав 


frequency polygons in figures 3(a) and 3(b). Ајвс shown in the figures 
are the normal distributions fitted in accordance with standard sta- 


tistical procedures by means of the calculated standard deviations of 
table I. | | 


The fit of the normal distribution to the data for airplane А 
shows some tendency of the normal distribution to underestimate the con- 
centration of values about the mean of the distribution while the fit , 
of the data of airplane В appears generally good. Іп order to test the | 
hypothesis that the observed samples are from a normal distribution, ! 
statistical procedures were applied to the observeü statistical param- 
eters. The results of these tests indicated that the distribution of 
An for airplane B could be considered from a normal distribution, whereas 
‘the distribution for airplane A would not be assumed from a normal dis- 
tribution. The magnitude of the departure from normality did not. nhow- 
ever, appear large in this case. | 


In order to permit examination of the behavior of the distributions 
at the larger values of load factor which are of concern, the observed 
relative cumulative frequency distributions and the fitted cumulative 
probability distributions are shown in figure 4. The curves for both 
the probability of exceeding given values and the probability of being 
less than given values are shown on semilogarithmic paper in order to 
permit comparison at both the large positive and negative acceleration 
increments. Examination of these figures indicates that the data for 
airplane В аге in excellent agreement with the fitted curves. For air- 
plane A, the over-all agreement between observed data and fitted curve 
appears reasonable although some discrepancy between the observed dis- 
tribution and the fitted curve is apparent particularly at the larger 
negative values of An. These discrepancies are in general, not large, 
however, and might, for example, be due to piloting-technique effects. 
However, the lack of consistency in the results indicates the need for 


further study of the question of the normality of the distribution of | 
load. : 
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In order to examine the question of normality somewhat further at 
this time, recourse was made to additional data on frequency distribu- 
tions of load increment under conditions similar to those represented 
in figures 3 and 4, Examination of these additional data indicates 
that although the distributions were generally close to normal, іп a 
number of cases, more large accelerations were experienced than expected 
for а normal distribution. Consideration of the time-history records 
from which these distributions were obtained indicated that these 
departures from normality were associated with lack of homogeneity in 
the turbulence intensity during the flight run. Only for short runs of 
constant level turbulence did the distributions appear normal. 


The manner in which the departures from normality arise in nonhomo- 
geneous turbulence may be clarified by ап example. Consider а flight 
run through rough air consisting of somewhat more severe turbulent con- 
ditions during the second part of the run. The over-all distribution 
of loads may for this condition be expected to consist of two normal 
distributions, one for each part of the run. Тһе two normal distribu- 
tions would be expected to have mean values of zero for the load factor 
increment but different standard deviations. The combined distribution 
for the whole run can be shown (appendix) to depart from normality with 
an excess of observations at the center and at the larger values of load 
increment. This mechanism can account for the departures from normality 
observed in the data examined. | 


In view of the foregoing indications, it would appear’ reasonable 
for the present to assume that for locally homogeneous turbulence of 
the type represented in figures 3 and 4 the frequency distribution of 
load increments may be expected to approximate a normal distribution. 
The approximate normality uf the load distribution for homogeneous | 
turbulence permits the use of the relations derived earlier between the 
frequency distribution of loads and the power spectrums and permits the 
representation of the load intensity by а single number, the standard 
deviation of the frequency distribution of load. This unification of 
the power-spectral-density function and the probability distribution of 
loads even for limited conditions is of considerable importance since in 
шапу load studies the relative loads in continuous rough air of two air- 
planes (or one airplane at two flight conditions) are of interest. For 
this case the use of the standard deviation as derived from the power- 
spectral analysis appears to provide a direct measure of the relative 
loads. . 


In view of the foregoing indications that a lack of turbulence homo- 


geneity can cause significant departures from a normal distribution of 
loads, further studies of the distribution of loads under various atmos- 
pheric conditions аге needed. It may be expected that the turbulence 
connected with such dynamic phenomena as thunderstorms which have a 


large-scale physical structure and short life cycle may not be stationary 
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random. For these conditions, the airplane load histories might have 

Е distributions that depart significantly from normality. This limita- 
tion may not, however, be too serious since even in these cases the 
standard deviation of loads may still provide a measure of load intensity 
although one not as simply interpretable in terms of the probability 
distribution as in the case of the normal distribution. 


Relation of Step Response Characteristics to Loads 
in Continuous Rough Air 


Examination of calculated responses to entry of a sharp-edge gust 
for airplanes representative of conventional transport types has indi- 
cated that the shape of the response curve up to the peak acceleration 
depends primarily on the gust penetration function (Küssner function). 
After the peak value which occurs close to 6 chords of gust penetration, ` 
the character of the response appears to be primarily a function of air-  - E | 
plane stability апа to approximate the short-period oscillation of the p 
airplane, Оп the basis of these properties, 8 limited series of- response . 
curves to a sharp-edge gust were selected to represent variations in о. 
airplane dynamic longitudinal stability. The incremental acceleration | 
responses in g's to а l-foot-per-second sharp-edge gust were assumed to 
: be given by the following expressions: 


LUE FIT m Ж-- - 


Аоба) = + sin db в 2222 05856) 
(31) 2 
1 - ы 
An(s) = = е7008-6)сов %.(в8-Ө (65852) 


Equations (31) represent а quarter sine wave up to а fixed peak value 


of 0.033g at 6 chords and a damped oscillatory function for the remaining 
portion of the response. The parameters b and kg сап be considered 


to represent the damping and frequency parameters, respectively, of the 

| response following the peak load. Іп figure 5, plots of equations (31) | 
are shown for nine selected cases covering three values for each of the 

two parameters. The values were selected to sample a wide range of'sir- 

plane response characteristics and represent variations in wave length 


from hO to 150 chords and variations in demping from light to almost 
critical damping. 
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In view of the linearity of the systems being considered, the 
multiplication of the response An(s) of equations (31) by a constant з 
results in the multiplication of the frequency-response function and 
the standard deviation of the output in turn by the same constant. Thus, 
the effect on the standard deviation of load of variations in the peak 
values of An for given values of Ко and b is apparent and the соп- 
29 sideration of an amplitude factor in equations (31) does not warrant 
| more detailed discussion. | 
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The frequency response function Т(іу) was determined analytically 


from the sharp-edge-gust response given by equations (31) by using equa- + 
tions (21). Тһе amplitude squared of' the transfer function | T(iv)|? is | 
given by i | | | : 


Е | (ау) [2 —— 


2 | 
(30) E 4 
У е; | 
(2 + Ко +у : | е, 
і 
| 
| | 
TN 
l оке Мају Фата 6v X 
(30)* i i | 
| 
E 3 | 
| j 
4 | 
сов - 4 | 
| у 12 j | 
| 
Plots of equation (32) are shown in figure 6 for the nine cases соп- : | 
sidered in figure 5. j 
| 
“4 
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The reduced power-spectral-density function of atmospheric turbu- 
lence, shown in figure 2, was used as the input. For any particular 
weather condition, the power-spectral-density function would have tc be 
adjusted by taking into accomt the actual average power of the input. 
This change would, however, only involve an appropriate linear change in 
the ordinate scale and would not affect the relative indications of the 
present comparisons. | | 
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Examination of figure 6 indicates that, ав the-damping of the 
oscillation is decreased (b decreased), the amplitude increases rapidly 
in the neighborhood of the oscillation frequency k,. The frequency 
response function remains relatively unchanged over the rest of the 
frequency range. On the other hand, variations in the oscillation fre- 
quency Ко, for fixed values of b, have a minor effect on the shape 
of the function but the peak values cHange (note changes in the ordinate 
scale) and occur at values of v close to ko. 
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The power-spectral-density function of the acceleration-increment’ 
output was obtained from equations (18) апа (22) by the relation 


_ ĝa) = 9,(0)| аса)|2 (33) | 


-- 


For convenience, ап airplane chord of 9.67 feet, the ваше value ав the 
chord of the airplane used to obtain the input spectrum, was assumed 
for this series of calculations. The power-spectral densities of the | 


acceleration-increment output S (9) for the conditions being con- i 
sidered are shown in figure T. Since the power-spectral-density func- 4 
tion of atmospheric turbulence was not known for frequencies less 

than 0.0016 radian per foot, the output spectrum could not be determined | 
in this region. The output spectrums were, however, extrapolated to ' 
zero at Я = О in order to complete the output spectrum at the low fre- 
quencies. The extrapolations are indicated in figure 7. 


С Examination of figure 7 indicates that as damping is decreased the 
' . output spectrums for each value of kg increase rapidly in the neighbor- 


hood of the oscillatory frequency @ = Ко = sar For given values of 
damping, however, the power spectrums do not vary appreciably in shape 
but the peak values shift in frequency with Ко. Thus, the total power 
of the spectrum as measured by the integral of the spectral function 

: appears to be largely independent of Ко and primarily а function of Ъ 
for the conditions investigated. This result can be seen from figures 6 
and 7 to be a consequence of both the variations in response functions 


. • and the rapid decrease of power of the gust spectrum with increasing 
и | frequency. 
| 
| | up — = - „== 
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The rapid increase of the loads with decreased damping is more 
clearly illustrated in figure 8, which shows the standard deviation of 
loads с (the square root of the integral of the power spectrum) as a 
function of frequency for the several values of b considered. The 
standard deviations of the acceleration increment 0 were obtained · 
from the output spectrums by using equation (29). Та the evaluation 
of equation (29), the output spectrums as extrapolated from a fre- 
quency of 0.0016 radian per foot to О (fig. 7) was used and the 
area of the output spectrums for values of Я greater than 0.028 was 
assumed negligible. The figure indicates clearly that the load level, 
as measured by the standard deviation, is largely independent of the 
frequency parameter Ко but varies appreciably when the damping param-. 
eter b is changed, the variation of the standard deviation being of 
the order of two to one for the conditions considered. 
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The results imply that the short-period response characteristics : 
and particularly the damping characteristics have an appreciable effect. | 
оп the airplane loads in continuous rough air. In contrast, the air- 
plane peak-load response to discrete gusts can be shown for the present 
example to be affected only to a minor extent by the short-period damping 
characteristics. In view of these differences between the indications of 
power-spectral and discrete-gust calculations, discrete-gust calculations ° 
may not adequately reflect differences in the gust loads in continuous 
rough air between airplanes differing in dynamic stability characteristics, 
In particular these results indicate that discrete-gust calculations may : 
not be adequate for the determination of loads in continuous rough air for i 
modern high-speed airplanes which in contrast to the airplanes of the past 
have relatively poor short-period damping. More complete accounting of 
airplane short-period characteristics in gust-load analysis would appear 
to be required in these cases. 
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The variations of the response to а step gust assumed in the present | 
illustration represent idealized conditions in which the characteristics 
of the airplane response to a step gust were changed in а simpJe manner. | 
in practice, the change of almost any airplane parameter will modify the 
response to a step gust іп в complex manner. Thus, a change of an air- 
plane parameter will affect the peak-load value, the location of peak, 
and both the frequency and damping of the subsequent oscillation. Con- 
sequently, the problem of optimum design for gust-load reduction is 5 і 
extremely complicated and beyond the scope of the present study. The | | 
effects of variations in airplane geometry on gust loads are, to some 
extent, indicated in the second illustration in which the complete | 
changes in the response to a step gust for limited variations in each 
of selected parameters of an idealized airplane are considered. | 
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· Effect on Loads of Some Variations in Airplane 


Geometric and Aerodynamic Parameters 


P 
анай 


In order to obtain an indication of the relation between some аіг- 
plane geometric and aerodynamic parameters and loads in continuous rough 
air, the power-spectral-density functions of load based on the turbulence 
spectrum of figure 2 were calculated for selected variations in airplane 
characteristics. An idealized transport airplane was used as a basic 
configuration. Characteristics of this basic configuration are listed 
in table II. Іп addition, values of each of 11 airplane parameters |2 
listed were varied separately in order to obtain an "increased" condi- ` 
tion and a "decreased" condition, as indicated in the table, (Їп the 
case of the center-of-gravity position, the increased condition repre- 
sents the forward center-of-gravity position.) The 23 conditions covered ) 
in these calculations represented stable airplanes with static margins 
varying from -0.252 to -0.506 where 68 is mean aerodynamic chord. 
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The responses to а unit step gust with а velocity of 1 foot per 
second A(x) for the two-degrees-of-freedom case, vertical motion and | 
pitch, were available for the 23 cases of table II from an unpublished 
Study based on the methods of reference 6. The frequency-response func- 
tions were determined by means of equations (17). Тһе power-spectral- 

: density function for acceleration output corresponding to the turbulence | 
input of figure 2 was obtained for each condition by use of equation (18). 

The standard deviation of acceleration increment was determined for each | 

| 

| 


сиве from equation (29) by the relation 


o 1 
С = П ом 


in performing this integration, it was assumed that the output for fre- U 
quencies of 2 greater than 0.028 was negligible. The output spectrums 

were also faired to a value of zero at 0 equal to zero as in the 

previous application. - 


/2 


As.a basis for comparing the results of the indications of the 
spectral analysis with the indications of conventional types of analysis 
based on responses to single representative gusts, a triangular gust. 
with a gradient distance (distance from zero to peak value) of 10 chords 
was selected as a representative gust condition. This gust condition 
represents an average gradient distance for the more severe gust loads 
and is frequently used for analysis purposes as a measure of the airplane 
loads in rough air. The реак load-factor value Anpay for a l-foot- | 
per-second triangular gust having a gradient distance of 10 chords was . 
calculated for the two-degrees-of-freedom case for each of the 23 condi- 
e tions considered. 
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. The results of the power-spectral and single-gust calculations are 
compared in figure 9 for each of the 23 conditions. The abscissa in the % 3 
figure represents the single-gust peak response, the ordinate the power- ` | | 
spectral-density measure of load intensity. For ease in comparing the i 
relative changes, both single-gust and spectral-calculation results are 
shown in terms of the basic airplane response. Thus, the abscissa is 
given as Armax /(A"max)pasic and the ordinate, as 9) Фав1с. The line 


Аппах С 0 
(Аллах) basic basic 


methods of analysis gave the same percentage change in load for a given 


change from the basic airplane, the plotted point would fall on the line 
of equality. 


~ а M ~ 
oa 


of equality 


is shown as a reference. If both 


Examination of the results shown in figure 9 indicates that to a j 
first approximation both the standard deviation of loads and the maxi- 
mum load for а 10-chord-gradient triangular gust show the ваше trends 
for variations of the airplane parameters considered. The largest 
changes in load level are associated with changes in airplane weight, 
wing area, air density, and slope of the wing lift curve as might be 
expected from the sharp-edge-gust relation. Inspection of the figure 
indicates that complex second-order differences that may be important 
exist in the indications given by the two measures of loads. Ав an 
example, for the variations of loads for changes of wing area the figure 
indicates that for the 20-percent increase in wing area (represented by 
the square), the value of шах | (Annex a increases about 18 per- 


cent. Тһе relative standard deviation o/dpasic, however, increases 


only about 7 percent. This difference appears to be а consequence of 

the greater sensitivity of the power-spectral analysis to the changes 

in airplane stability introduced by the increased wing area; the increase 
in wing area results in а small decrease in the frequency of the response 
to a step gust but a more pronounced increase in the damping. The present 
results also suggest similar differences between the two analyses for 
changes in airplane weight and slope of the wing lift curve. 


j 
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The variations considered for the remaining parameters such as tail 
length, tail area, and tail slope of the lift curve, in general, yield 
minor variations in loads. The variations in ^na x/ (Anmax baigie for 


these cases are however less than #2 percent while the variations for the 
standard deviations are generally somewhat larger, il percent. In fig- 
ure 9 the l2-percent-chord rearward movement of the center-of-gravity | 
position ly yields about а 3-percent increase in Anmax |(Аопах — 


but а 9 percent increase in the relative standard deviation. In this 
cage, the larger increase in the standard deviation appears largely 
associated with the movement of the peak of the frequency response func- 
tion to lower frequency and thus larger gust spectral power. 
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А spectral analysis and the single~gust реак-1сай response yield trends 
- for changes іп parameter values that are generally consistent. Differ- 
ences that mey be important are, however, discernible and appear associ- 
ated with the greater dependence of the standard deviation of loads on 
the dynamic stability characteristics of the airplane. It would be 


expected from the first application that for larger variations in the | 
{ 


For the stable airplane configurations considered, the рочег- | 
4 
4 


damping characteristics of the short-period oscillation, the differences 
between the two types of analysis would be appreciably larger than 
obtained in the present illustration. 


" SUMMARY OF RESULTS 


The anelysis of the application of power-spectral methods of 
analysis to gust-load problems has indicated the following results: 


1. The application of power-spectral methods of analysis to load - | 
calculations provides a measure of load intensity for continuous rough 
- air in terms of the standard deviation (root mean square) of the load ! 


Е output • 


5 | T 2. The probability distribution of load intensity in homogeneous 
rough air appears to approximate a normal distribution. | 


3. Рог the case of the normally distributed output, the standard 
deviation of load completely describes the probability distribution of 
loads specifying the proportion of total time that various load values 
are exceeded. | 


4 “К-ға м --- 
a 
- 
= 


| 
1 
| 
| 
4, The application of power-spectral relations to a selected series 2 ж 
of systematic variations in the frequency and damping characteristics 
of the airplane load response to a step-gust input indicates that the 
damping characteristics of the indiclal response are the primary charac- | 
teristics in determining the loads in continuous rough air. This result : 
appears of significance in regard to high-speed airplanes and missiles 
where the short-period longitudinal damping may be poor. | 
5. Calculations for a limited series of conventional and stable 
airplane configurations indicates that the loads in continuous rough 
air for variations in individual airplane geometric and aerodynamic 
parameters are to a first approximation adequately reflected by the 
| peak-load response to the arbitrary 10-chord triangular gust commonly 
E vned. However, differences are discernible between the indications of 
5 the continuous-gust and discrete-gust calculations and appear largely 
"a associated with differences in the effects of stability changes on the 
: loads for the two gust conditions. 
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SUGGESTIONS FOR FUTURE RESEARCH 


| | 

: 22 Тһе foregoing results appear to indicate that power-spectral 
— methods are well-suited to the calculation of loads in continuous rough 1 
air. A number of problems exist end require further investigatione | | 
These include | | 
Li 


(1) The determination of the power-spectral-density functions of 
atmospheric turbulence for a wide range of atmospheric conditlons. 


(2) The determination of the conditions under which normal distri- | 
butions of the load time history apply. 


(3) The investigation of the effects of the nonuniformity of gusts 
across the airplane span on the application of these techniques. 


(4) The experimental verification of theoretical determinations of | 
airplane indicial responses over the range Scared to determine usable ' E 
frequency-response functions. 


7 (5) The investigation of the relations between-the power крек йа 
and the frequency distribution of load-increment Peake and other quanti- 
ties of interest in structural design. . . 


- 


qe E JM"ECcU Mo. 


—-— 


=з 


Langley Aeronautical Тавата гу, 
National Advisory Committee ПЕРЕТЕРТИ 
Langley Field, Va., September 29, 1952. 5 
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| APPENDIX | 
DISTRIBUTION OF GUST LOADS IN NONHOMOGENEOUS ROUGH AIR 


If the two normal frequency distributions of load increment with 
zero mean are defined by 


Wet pe 


a 
Pe ee ee и ee ВУ ДК јави 
e 


N,f (An) = ~ e 
o Non А 
(А1) d 
: Mr | 
No 210 
Nf (An) = e Е | 
; - | a5 Von | 
| 
у where i | | ‚- 
21 (Аа) ,f5(4n) respective probability density functions i | 
№1 ‚№2 respective number of observations for each 
distribution Ab 
0) 95 standard deviations of respective distributions . ! 
| 
F the combined distribution is defined by | | | 
| E 
(A2) 


Nf(An) = №121 (Ап) + Nofo(dn) 


where N = Nj + No. This distribution is examined for normality. 


Since the two frequency distributions Njfj(^n) and Nofo(An) аге 


each normal with mean zero, the combined distribution Hf(^n) is also 
symmetrical. about zero. Consequently, for the moments of the frequency 
- "distribution ы, defined by : 


* Ds 


~ со 


where m is a positive number and designates the moment order, the 
following relations apply 
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Ва = | ” (да) (аш) alan) | (АЗ) NE | 
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ug Ё О (m even) 


For а normal distribution, the coefficient of kurtosis a, defined 


by 
J | 
u H 7 T ES | 
I = cn = + (A5) | | 
Ho c | | | | 
- - X 
has а fixed value | 
а, = 3 · (Аб) | 
, | 


In order to examine the normality of the distribution defined by 
equation (A2), consider only the value of oy, the coefficient of 
kurtosis, for that frequency distribution. 


By definition, the coefficient of kurtosis a, for the distribu- 
tion of equation (A2) is given by | 


= |" (Ал) 17) а(Ап) + (| i (да) “N fo (Да) а(Ап) 


а 
| A (An)? hr e (An) + Nota (^) | (An) 


(А7) 


Ny + № 


3 l 4\ 
Ny + No (А + №02 ) 
e = 


E (48) 


2 оҳ: 
№01“ + н202 


Ny + No 


From equation (A6), a necessary condition for the distribu- 4 
tion Nf(An) to be normal is that а, be-equal to 3. From equa- | | 


tion (A8) this condition is obviously true only when | 


{ | 
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Which when simplified yields | т | 
| 
4 
| 
х | 
| 


в. p 1 2 2\2 | 
E "uq uo) E (ве; + Ноор (А9) | ) 
ы Expanding the term on the right and simplifying yields = Et 
i | 
| E 
Equation (A10) is thus а necessary condition for the distribution А 
defined by (А2) to be normal, However, from the inequality 1 | 
| 
a^ + b^ > гађ „(а ў v) (A11) j 
| | | 
| the only condition for which equation (A10) holds is.when су = dp. l i 
; Thus the distribution f(An) has а normal distribution only for the i 


trivial case бу = 02. If су ў ор, it also follows from equation (A10) 
that the coefficient of kurtosis a, for Nf(An) given by equation (A8) 
is greater than 3. The combined distribution consequently has an excess 
of kurtosis. Before considering the significance of the excess kurtosis, 
it 1s well to note that the present derivation was restricted to the 
combination of two distributions for simplicity. It ів simple. although 


* tedious to verify that similar results are obtained for the combination 
of three or more distributions. i 
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The significance of the excess kurtosis for the combincd distribu- 
tion may perhaps be best indicated vith a sketch. The sketch shows two 
symmetrical relative frequency dis- 
tributions having the same standard 
deviations, one a normal distribu- 
tion (oy = 3) the other a nonnormal 
distribution with a, > 3. The non- 
normal distribution is derived from 
the normal by а shift of mass out- 
ward from the central portion of the 
distribution and also a shift of 
mass inward to maintain the same 
standard deviation for the two савев, 
Thus, it is apparent that the pre- 
dominent feature of the distribu- 
tion having excess kurtosis is a 
greater concentration of mass at the 
center and at the outboard regions 
of the distribution than for the · 
normal-distribution case. 
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TABLE I.- FREQUENCY DISTRIBUTIONS OF ACCELERATION INCREMENT 


Airplane A Airplane B | 


Acceleration increment, Acceleration increment, 


-0 2177 to 0.727 


У У Nec m заваан - 
. 
. 


ч, 


-. 124 to 
-.6TT to 
~,627 to 
-.577 to 
-. 521 to 
-.477 to 
~-.427 to 
-.3TT to 
-. 321 to 
-.27( to 
-.c27 to 
-.177 to 
-.127 to 
-.077 to 
-.027 to 
~023 to 
«073 to 
‚123 to 
• 173 to 
‚223 to 
273 to 
„323 to 
•373 to 
423 to 
+73 to 
. 223 to 
.• 273 to 


~. 677 
-.627 
-.5TT 
-. 527 
-. 571 
- A427 


--.3Т71 


-.327 
-. 211 
-. 221 
-.177 
-.127 
-.077 
-.027 
2023 
• 073 
• 123 
• 1713 
•223 
„2{3 
„323 
• 373 
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TABLE II.- VALUES OF AiRPIANE PARAMETERS 


fu = 1 ft/sec; V = 308 ft/sec; cy = 7.58 rt] 


Sy, per radian . . « о. 
at, per radian 
S, sq.ft... 
5+, ва ft e о 


р, slugs/cu ft 


"" “Values indicate parameter changed Рог given condition 
whereas other parameters remain same as basic condition. 


Center of gravity rearward. 
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pheric 


per second, 


Figure 1.- Normalized power-spectral-density function for atmos 
vertical gust velocity Фј(ш) for airspeed of 300 feet 
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-— o Figure 2.- Normalized power-spectral-density function for atmospheric 
| vertical gust velocity $, (R). 
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(a) Airplane А. 
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Figure 3.- Comparison of observed frequency distribution and fitted normal 
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(b) Airplane B. 
Figure 3.- Concluded. 
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(а) Airplane A. 


Figure №.- Comparison of observed relative cumlative frequency distribution 


with fitted normal probability distribution. 
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(b) Airplane B. 


Figure №.- Concluded. 
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edge gust for several stability conditions. 
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(b) kg = 0.07. 
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(c) Ко = 0.10. 
- Modulus squared of transfer functions for selected variations 


the frequency parameter Ко and the damping parameter b. 
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Figure 6. 
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(a) Ко = 0.0}. 
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Figure 7.- Power-spectral-density function of load-factor output. 
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Figure 8.- Standard deviation of loads as a function of the frequency 
parameter К, for various values of the damping parameter b. 
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Figure 9.- Comparison of standard deviation of load 


with peak load 
Flagged symbols 


-gradient triangular gust. 


denote below-basic parameters. 


response to a 10-chord 
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